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In previous work we showed that bright atomic solitons 
can arise in spinor Bose-Einstein condensates in the form of 
gap solitons even for repulsive many-body interactions. Here 
we further explore the properties of atomic gap solitons and 
show that their internal structure can be used to both excite 
them and control their center-of-mass motion using applied 
laser and magnetic fields. As an illustration we demonstrate 
a nonlinear atom-optical Mach-Zehnder interferometer based 
on gap solitons. 



I. INTRODUCTION 

The experimental demonstration of Bose-Einstein con- 
densation in atomic vapors has rapidly lead to spectac- 
ular new advances in atom optics. In particular, it has 
enabled its extension from the linear to the nonlinear 
regime [Q, very much like the laser led to the develop- 
ment on nonlinear optics in the 1960s. It is now well es- 
tablished that two-body collisions play for matter waves 
a role analogous to that of a Kerr nonlinear crystal in op- 
tics. The first experimental verification of this analogy 
was four- wave mixing In addition, it is also possible 
to nonlinearly mix optical and matter-waves, as demon- 
strated recently in matter-wave superradiance and 
in the first realization of a phase-coherent matter-wave 
amplifier In addition, it is known that the nonlin- 
ear Schrodinger equation which describes the condensate 
in the Hartree approximation support soliton solutions 
1^. For the case of repulsive interactions normally en- 
countered in BEC experiments, the simplest solutions 
are dark solitons, that is, "dips" in the density profile 
of the condensate [f7|-pO{ . These dark solitons have been 
recently demonstrated in two experiments [|ll|,^ which 
appear to be in good agreement with the predictions of 
the Gross-Pitaevskii equation. 

While very interesting from a fundamental physics 
point-of-view, dark solitons would appear to be of lim- 
ited interest for applications such as atom interferometry, 
since what would be desirable is to achieve the dispersion- 
less transport of a spatially localized ensemble of atoms, 
rather than a "hole." In that case, bright solitons are 
much more interesting. However, the problem is that 
large condensates are necessarily associated with repul- 
sive interactions, for which bright solitons might appear 
to be impossible since the nonlinearity cannot compen- 



sate for the kinetic energy part (diffraction) in the atomic 
dynamics. While this is true for atoms in free-space, this 
is not the case for atoms in suitable potentials, eg. op- 
tical lattices. This is because in that case, it is possi- 
ble to tailor the dispersion relation of the atoms in such 
a way that their effective mass becomes negative. For 
such negative masses, a repulsive interaction is precisely 
what is required to achieve soliton solutions. This re- 
sult is known from nonlinear optics, where such soliton 
solutions, call ed gap solitons, have been predicted and 
demonstrated |13 IJ]. 

We showed in previous work that such bright soli- 
tons are also possible in matter-wave optics using spinor 
condensates, but did not discuss explicitly how to excite 
and control them. The present paper addresses these 
questions, and shows that a combination of optical and 
magnetic fields can be used to generate solitons of var- 
ious velocities, and subsequently, to control them, split 
and recombine them, etc. This magneto-optical control 
results from the use of spinor condensates, and can be 
achieved at minimal cost in terms of atomic loss. In 
particular, we illustrate how to realize a Mach-Zehnder 
interferometer for bright atomic solitons. This opens up 
the way to intriguing new ways to manipulate and trans- 
port coherent matter in ways which complement those 
offered by standard optical tweezers or by atomic wave 
guides (l|0. 

Amongst Dan Walls many contributions to nonlinear 
and atom optics he was one of the first to recognize the 
emerging area of nonlinear atom optics, and he intro- 
duced the idea of atomic solitons traveling in laser beams 
already in 1994 along with Zhang and Sanders ||l^. The 
current wave of interest in atomic solitons therefore finds 
its origin in Dan's seminal work, and we all greatly miss 
the originality, insight, and energy he brought to the field. 

The paper is organized as follows: Section || briefly 
reviews the analysis leading to the predictions of gap 
solitons [fsi, while section H discusses some of their 



most important characteristics. Section IV exploits these 



properties to develop tools to excite gap solitons and con- 
trol their dynamical behavior, leading to the demonstra- 
tion of an atomic Mach-Zehnder interferometer in Sec- 
tion Finally, section VI is a summary and outlook. 



While the main text uses a simple model of light-matter 
interaction for notational clarity, a more realistic cou- 
pling scheme, using the full hyperfine structure of the 
Sodium 3S'i/2-3P3/2 transition, is presented in Appendix 
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II. PHYSICAL MODEL 

To set the stage for our analysis, we first briefly review 
the main ingredients of the theory of gap solitons of Ref. 

. The system we consider consists of a Bose-Einstein 
condensate interacting with two counterpropagating, fo- 
cused Gaussian laser beams of equal frequency ui but 
opposite circular polarizations, see Fig. |l|. The optical 
dipole potential associated with the applied laser beams 
is assumed to provide tight transverse confinement for 
the BEG in the {X, Y) plane, thereby forming a cigar 
shaped condensate of transverse cross-sectional area At- 
In the following, we confine our discussion to the one- 
dimensional dynamics of the BEG along the Z-axis for 
simplicity. 

In addition to supplying a transverse optical potential, 
the laser beams can drive two-photon transitions between 
different Zeeman sublevels of the atomic ground state. 
For illustrative purposes we consider the case of Sodium 
and the two-photon coupling of the Zeeman sublevels | — 
1) = \Fg = l,Mg = -1) and |1) = \Fg = l,Mg = 1). For 
example, starting in the | — 1) state this process involves 
the absorption of a (T+ photon from the right propagating 
laser beam followed by emission of a (t_ photon into the 
left propagating laser beam. 

We must of course assume that the excited states in- 
volved in the atom-field interaction are far-detuned from 
the applied laser frequency, a necessary requirement to 
avoid the detrimental effects of spontaneous emission. 
This however raises a serious issue, since for the alkali 
atoms that we have in mind the two-photon coupling 
strength vanishes in the limit when the detuning is large 
compared to the excited-state hyperfine splitting ||l9|,^ . 
But this difficulty can be circumvented: We show in the 
Appendix that by using a four-photon scheme involving 
an additional 7r-polarized field incident perpendicular to 
the Z-axis, one can achieve an effective two-photon cou- 
pling between the states | ± 1) which survives in the limit 
of large detunings. We therefore proceed with our two- 
photon coupling model and refer the interested reader to 
the Appendix for details. 

By restricting our attention to the coupled states | ± 1) 
the effective single-particle Hamiltonian for our model 
system can be written as [|l5|,pi[ 



Heff = 



p2 
f_Z_ 

2m 



ghS'[\l) 



„2iKiZ 
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where we have omitted constant light-shift terms. Here 
Pz is the atomic center-of-mass momentum operator 
along Z, m the mass of the atom, Ki = / c is the magni- 
tude of the field wave vector along Z, g is a coupling con- 
stant between the ground and excited states characteris- 
tic of the atom and transition involved, S' = V^S^ /h^S, 
with the detuning S = uJi — uja, £ is the laser field ampli- 
tude, and V is the reduced dipole moment for the 35i/2~ 
3P3/2 transition. 



The first term of the effective Hamiltonian de- 
scribes the quantized atomic center-of-mass motion, and 
the remaining terms give the effective coupling of the two 
Zeeman sublevels via the applied laser fields. The expo- 
nential terms exp(±2iif;Z) arise from the fact that the 
two-photon transitions involve the absorption of a photon 
from one light field and reemission into the other. Finally, 
introducing a spinor macroscopic condensate wave func- 
tion 1'{Z,t) = [i'i{Z,t),^_i{Z,t)]'^ normalized to the 
number of atoms N, and including the many-body ef- 
fects via a mean-field nonlinearity, we obtain the coupled 
Gross-Pitaevskii equations 



(2) 



asc is the s-wave scattering 
^, and we have assumed that 



where U = A-kTl ascl°fn^ 
length, l^-p = 1^-112 + 
the magnitude of the self- and cross-nonlinearities are 
equal for simplicity. 

It is convenient to re-express Eq. (H) in dimensionless 
form by introducing the scaled variables r — t/tc, z — 
Z/lc and ijjj — jl \/p^ where 



_i_ _ UfiKi 
go' m 



ghS' 



U 



(3) 



Equations (0) then become 



sgnigS'/U)\^\'(^^^ 



■01 



(4) 



where M ~ gS''m/2hKf is a mass-related parameter such 
that ki = Kilc ~ \/2M. Throughout this paper we use 
a characteristic laser intensity of / = 50 W/cm^ for each 
of the two counter-propagating laser beams and a wave- 
length of A/ = 985 nm. In addition we use g = —1/4, 
which is simply the Clebsch-Gordan coefficient for the 
Fg = 1 ^ Fe = 1 transition. This results in the charac- 
teristic scale values tc = 685 ^s, 
8.59 HTcr^ and M = 0.0065. 



Ic = 12.1 ^m, pc = 



III. GAP SOLITONS 

The spatially modulated coupling between the optical 
fields and the condensate induces a single-particle band 
structure with regions of negative effective mass. As men- 
tioned in the introduction, this leads to the possibility 
of bright atomic solitons even for repulsive interactions 
p^ . Their energy lies in forbidden gaps of the linear 
band structure, hence the name gap solitons. 
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A. Analytical Solutions 

Approximate analytic expressions for the gap solitons 
can be obtained by expressing the spinor condensate com- 
ponents as 



extract further important soHton properties, such as the 
number Ns of atoms in the sohton shown in Fig. H, and 
the sohton width Ws = Wslc shown in Fig. ||. Specifi- 
cally, the number of atoms in a particular gap soliton is 
given by 



^±ikiz^-iT/iM 



0±i(2,t), 



(5) 



where the field envelopes (j)±i are assumed to be slowly 
varying in space compared to l/fcj. Neglecting then the 
second-order spatial derivatives yields the coupled partial 
differential equations 



dr dz 



1 

1 



(6) 



where the choice ±1 — sgn{gS' /U). For a red-detuned 
laser and our choice of g this becomes ±1 = sgn(C/). 
Aceves and Wabnitz |22| have shown that these dimen- 
sionless equations have the explicit two-parameter gap 
soliton solutions (see also Ref. |14]) 



= ± 



sin(7y) 
(3 sin(ry) 



„2B 



g2e _|_ g±i,7 



sech 



2 
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c2e 



_|_ g±ir, 



sech 9 ± 



irj 



with — l<u<lisa parameter which controls the soliton 
velocity, < 77 < tt is a shape parameter, and 



1 



1 



(8) 



= — 7sin(?7)(z — wr), cr = —7 cos(?7)(wz — t). (9) 

Since we are interested in creating bright solitons in the 
presence of repulsive interactions we restrict ourselves to 
sgn(C/) = -1-1, corresponding to the choice of the upper 
sign in the analytic solutions. 

The characteristic length scale associated with the soli- 
tons is Ic, so that the approximate solitons (Q) are valid 
for Kilc = 1/2M >> 1. This inequality is well satisfied 
for our choice of parameters, which gives l/2Af = 76.9. 



B. Characteristic Properties 

From the dependence of the hyperbolic-secant on 9 ~ 
— 7sin(7])(z — vt) in Eqs. (^), we identify the gap soliton 
parameter v — Vg/Va as the group velocity Vg of the soli- 
ton in units of the recoil velocity Vr ~ Ic/tc — TiKi/m. 
Since — 1 < v < 1, the magnitude of the group velocity is 
bounded by the recoil velocity, which is Vr — 1.77 cm/s 
for the present case of Sodium. From Eqs. (Q), one can 



At / dZ[|vI/i(Z,i)|2 + |vI/_i(Z,t)|2], 



(10) 



(7) 



where At is the effective transverse area. Fig. ^ il- 
lustrates two general soliton properties, namely that in- 
creasing the shape parameter rj increases the atom num- 
ber, but that faster solitons have lower atom number. 
Similarly, Fig. |^ shows that increasing the shape param- 
eter 77 decreases the soliton width and that faster solitons 
have narrower widths. 

Given that the analytic gap soliton solutions (^ hold 
for broad envelopes, we confine our attention to the case 
ry < 1, with atom numbers in the range Ng ~ 10^ — 10^, 
and soliton widths Ws ~ 40 — 100 /im for the parameters 
at hand. 

We can gain further insight into the structure of the 
gap solitons by taking the more extreme limit 77 << 1 
of Eqs. (0). Using the definitions (||), and returning to 
dimensional units we have then 

/37 V 2 

^jJZ,0) = ^./^sec/i(Z/Wo)(-l)"e'''^"^'^^+"', (11) 
7 V 2 

where 



W, = 3AAWn = 3.44 



IcVT 



i^ = -f (12) 



Here Ws — 3.44Wo is the soliton width, the factor 3.44 
being the numerical conversion from the width of the 
hyperbolic secant to the 1/e^ width of the distribution, 
and K — k/lc & velocity-dependent wave vector shift. 
This expression agrees well with the features displayed 
in Fig. H in that the width decreases with increasing 77 
and V. The soliton atom number obtained by combining 
Eqs. (0) and ([ll]), 



Ns = 2{AtIcPc)v ■ (1 



(13) 



correctly predicts the scaling properties of Fig. ^ 

An essential point to keep in mind is that the gap soli- 
tons are coherent superpositions of the two Zeeman sub- 
levels, and the approximate solutions ( pi] ) contain im- 
portant information on the phase and amplitude rela- 
tions that need to be created between them to success- 
fully excite and manipulate gap solitons. In particular, 
they show that there is always a spatially homogeneous 
TT phase difference between the two states. In addition, 
the two components have the spatial wave vectors 



K±i = K±Ki. 



(14) 
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with K = —^v/lc- For our parameters, Ki = 6.38 /im"-*^ 
and \K\ < 1/lc — 0.08 /im~^. Despite the fact that it 
is so small, \K\ is an important factor since it controls 
the soliton velocity. Finally, it follows from dividing the 
amplitudes of the two components that 

^^ijZ^t) 1 ^ fl + v\ 

which shows that their relative occupation depends on 
the soliton velocity parameter v. For v = the sublevels 
are equally populated, but as u ^ 1 the |1 > sublevel 
has a larger population, and vice versa for v —^ —1. 

The characteristic time scale for the evolution of the 
gap solitons can be determined from the plane-wave ex- 
ponential factors in Eqs. (|^). Converting to dimensional 
form the soliton period ts is defined as the time to accu- 
mulate a 2tt phase, or in the limit 77 —s- 

ts = 2TTtc\/l-v'^. (16) 

Physically, corresponds to the internal time scale for 
the gap soliton. In order to observe a soliton-like be- 
havior, it is therefore necessary to investigate the atomic 
propagation over several periods. 

We conclude this section by noting that Aceves and 
Wabnitz have shown that the gap soliton solutions 
(^ are stable solutions of Eqs. (||) in that they remain 
intact during propagation, even when perturbed away 
from the exact solutions. However, one should remember 
Eqs. ^ are only an approximation to the exact system 
of Eqs. (^ so that in general the gap solitons are solitary 
wave solutions only. As such, they are not guaranteed to 
be absolutely stable. 

IV. GAP SOLITON CONTROL 
A. State manipulation 

Summarizing the previous section, gap solitons require 
the right population in each Zeeman sublevel, a phase 
difference of n between these sublevels, and appropriate 
plane- wave factors e*^±^'^. The shapes of the Hartree 
wave functions corresponding to the two Zeeman sub- 
levels are hyperbolic secant, which we approximate by 
a Gaussian in the following. They could for example 
be initialized in an optical dipole trap ||2^ . Manipulat- 
ing the gap solitons is therefore reduced to the problem 
of controlling the populations and phases throughout the 
spinor condensate. This is achieved via a magneto-optical 
control scheme involving a combination of pulsed coher- 
ent optical coupling and of phase-imprinting using spa- 
tially inhomogeneous magnetic fields. 

The coherent optical coupling can be achieved e.g. by 
a laser pulse of frequency toi propagating perpendicularly 
to the Z-axis and with linear polarization perpendicular 



to that axis. For sufficiently short pulses, one can ne- 
glect changes in the center-of-mass motion of the atoms 
during its duration, leading to a very simple description. 
We assume for simplicity a plane- wave rectangular pulse 
of duration tp and of spatial extent large compared to 
the soliton. The Hamiltonian describing the coupling 
between thispulse and the condensate is then the same 
as in Eq. (|l|) without the linear momentum exchange 
terms exp{zL2iKiZ) and the kinetic energy term, and 
with S' — > S'p. The state of the system after the pulse 
is then easily found to be 

f *i(^p) \ cosx ismx \ ( *i(0) \ 

1^ *-i(ip) ) \ -isinx cosx J \ *-i(0) J 

. M.,„ (*.;»)) ^ ,17, 

where x ~ is the excitation pulse area and the op- 

erator Mi can be used to control the population transfer 
by an appropriate choice of x- 

The required phase relation between the two states 
can be achieved via Zeeman splitting. Considering for 
concreteness a spatially inhomogeneous rectangular mag- 
netic field pulse of duration is we have, neglecting again 
all other effects, 

ih-^^^±i{Z,t) = ±^iBgF {Bo + B'Z) *±i(Z,t), (18) 

where (7f is the Lande g-factor of hyperfine ground state, 
/is is the Bohr magneton, Bq the spatially homogeneous 
component of the magnetic field, and B' its gradient, the 
direction of the magnetic field being along the Z-axis. 
The application of this field results in the state 

f ^litB) \ _ /e'(''+^-^) \ / ^i(O) \ 

M/_i(tB) J ^\ J [ vl/_i(0) J 

^Msi^,Ks)^^^lll^y (19) 

where 1) = —{^.B9F/K)BatB and Kb = —{^.Bgp /K)B'tB 
are the imprinted phase shift and phase gradient (or wave 
vector), respectively. That is, the application of the mag- 
netic pulse results in a phase difference of 2d between the 
two Zeeman sublevels, and in addition it imparts them 
wave vectors ii^B- 

We note that although we have treated the pulsed ex- 
citations above in an impulsive manner to illustrate their 
action, our simulations describe the actions of the pulses 
correctly to the equations of motion (^. The numerics 
confirm the accuracy of the impulsive approximation for 
the parameters at consequence of the fact that we con- 
sider pulse durations significantly shorter than the soliton 
period ( pj| ) of ts = 4.3 ms. 

B. Gap soliton excitation 

To illustrate how stationary and moving solitons can be 
excited using the proposed magneto-optical state scheme. 
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we start from a scalar condensate in the | — 1) state, 
*(Z,0) = [0, *o(-2)]^, with spatial mode 

with Ws and the width and atom number of the gap 
soliton desired, see Figs. || and ||. This Gaussian is cho- 
sen to approximate the hyperbolic-secant structure of the 
analytic gap soliton solution in Eq. For a station- 

ary solution we need to prepare the Zeeman sublevels 
with equal population and with a tt phase difference. We 
further need to impose wave vectors which are equal in 
magnitude but opposite in sign K±i = ±Ki. This can 
be achieved by applying a laser pulse with area x = ''''/4, 
followed by a magnetic pulse with ■& = tt/A and Kb = Ki. 
The state then transforms as {t ^ tp + ts) 

^{Z,t) = MBin/'i,Ki)ML{T:/A)^iZ,0) 

= ^(,-4-.K,z)*o(Z). (21) 

Using the atomic parameters of Sodium, this situa- 
tion can be realized for a 10 /^s light pulse with / = 
2.69 KW/cm^, and a 200 /is magnetic field pulse with 
Bo = 0.89 mG and B = 72.5 G/cm |6|. Fig. |shows the 
resulting stable evolution of the total density |\1/(Z, i)|^. 
As a result of the Gaussian approximation to the exact 
solution there are some slow oscillations imposed on the 
motion, but the solution remains centered at Z = and 
stationary over a time t = 200 ms, much longer than the 
soliton period ts ~ 4.3 ms. 

The excitation of a moving soliton is slightly more 
complicated since the velocity dependent wave vector 
K = —jv/lc in Eq. ( pT| ) is no longer zero. To deal with 
this situation, we first impart the wave vector K to the 
initial | — 1) state in Eq. ( [20[ ) using a magnetic pulse with 
i) = 0, Kb = —K, which for Sodium and a 200 mag- 
netic pulse can be realized using Bq = 0,B' = -1-0.094 
G/cm. Here we specifically take v — 0.1, so that 
K — —8306 m~^. We further recall that moving soli- 
tons must have unequal populations of the two Zeeman 
sublevels. For v = 0.1, we have from Eq. ( [l5| ) that 
the ratio between the |1) and | — 1) populations should 
be 1.22. This is achieved by a coherent optical coupling 
with X — 0.835, and / — 2.69 KW/cm , correspond- 
ing to a pulse duration of 10.6 fis. Finally, we impart 
the wave vectors ztKi and the tt phase difference be- 
tween the states the Zeeman sublevels with a magnetic 
pulse, as for the stationary soliton above. Summarizing, 
the full magneto-optical control sequence is described at 
{t = 2tB + tp) by 

■^{Z,t) - Mb{ti/A,Ki)Ml{x)Mb{Q.~K)-^{Z,Q) 

= ( cosxe' -e^(--<-) ) ^''^ 

Fig. 1^ shows the resulting numerical simulation of a gap 
soliton with v = 0.1 [Vg = 0.18 cm/s), illustrating its 



stable propagation is exhibited over many soliton periods. 
We have used the same scheme to launch gap solitons over 
the full range of velocities. 

C. Soliton splitting 

The next application of magneto-optical control that 
we consider is soliton splitting. To achieve this goal, we 
take advantage of the fact that for fast solitons almost all 
of the population is in one Zeeman sublevel, and the other 
state can be viewed as a small perturbation. For example, 
for V = 0.5 the ratio of the populations is already 6. This 
implies that the relative phase between the two states is 
no longer of importance, so we need only concentrate on 
getting the plane- wave factors right. 

Assume for concreteness that we start from an initial 
condition with Nq atoms in the | — 1) state, and apply 
a laser pulse of area 7r/4 to transfer half the population 
to the |-|-) state. We can then apply a magnetic pulse to 
impose the wave vectors ■iz{Ki -f A'), to the | ± 1) states, 
with K — —^v/lc corresponding to a given velocity \v\. 
Now if we choose Aq = 2 A^s , with A^^ the atom number 
for that velocity, then for |?;| large enough we may ex- 
pect to see oppositely moving solitons emerge from the 
initial state. This is illustrated in Fig. ^, which clearly 
illustrates the emergence of two solitons with opposite 
velocities. Fig. ^ which shows the density profiles for 
the individual Zeeman sublevels for a time t = S2 ms, 
confirms that that each gap soliton indeed comprises a 
dominant Zeeman sublevel plus a small component of 
the other state. During the early stages of propagation, 
the solitons rearrange their phase and shape before set- 
tling down. This is accompanied by some slowing down 
and the familiar shedding of "radiation." The emerging 
solitons are therefore slower than their "design velocity." 
For the specific example of Fig. ^, the actual velocity is 
found to be 1.08 cm/s. 

D. Soliton reversal 

In addition to offering the possibility of exciting mov- 
ing gap solitons, magneto-optical control can also reverse 
their direction of propagation along Z . This can again be 
achieved by making use of the fact that for fast solitons 
almost all the population is in one Zeeman sublevel, so 
that we need only concentrate on getting the plane- wave 
factors right. In the numerical simulation of Fig. ^ we 
create a gap soliton with v ~ 0.59 [Vg = 1 cm/s) and let 
it propagate for 100 ms. This soliton consists essentially 
of a plane-wave factor Q^i^+i^i)^ multiplying state |1). 
At t = 100 ms we apply a laser pulse of area of 7r/2 to 
transfer all the population to state | — 1) and change the 
phase to e~*(^+^')^ using a magnetic field pulse. This 
results in the same soliton as before the control sequence. 
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but with opposite velocity. Note however the loss of some 
atomic population to "radiation" in the process. 

V. ATOMIC MACH-ZEHNDER 
INTERFEROMETER 

As an illustration of the potential use of gap soli- 
tons employing magneto-optical control here we consider 
a nonlinear atomic Mach-Zehnder interferometer. Soli- 
tons present some advantages for atom interferometry in 
that they are many-atom wavepackets which are immune 
to the effects of spreading, hence allowing longer path 
lengths, and also reduced signal-to-noise for large atom 
numbers. Typically many-body effects limit the utility of 
high-density wavepackets due to spatially varying mean- 
field phase shifts, but solitons have the cardinal virtue 
that they have fixed spatial phase variations (for our sit- 
uation this applies for faster solitons). Thus solitons may 
provide a key to making maximal use of high density 
sources for atom interferometry. Indeed they have long 
been advocated for all-optical switching applications due 
to these very properties. 

Our specific demonstration of a nonlinear atomic in- 
terferometer based on solitons involves an initial scalar 
condensate that is split into two oppositely moving soli- 
tons along Z, see Fig. ^ for i < 60 ms. At i = 60 
ms laser and magnetic pulses are applied which act to 
reverse the direction of the two solitons. The process of 
reversing causes some loss of atoms in both solitons as be- 
fore. The two reversed soliton components come together 
again at i = 120 ms. Since the colliding solitons are pre- 
dominantly in opposite orthogonal Zeeman sublevels, the 
interference pattern appearing during the collision is due 
to the contamination of each soliton by the other state, 
see e.g. Fig. Fig. |l^ shows the interference at the 
soliton collision in the total density (solid line) and also 
the individual Zeeman sublevels, with a fringe contrast 
of around 30%. These results demonstrate the potential 
use of gap solitons for realizing nonlinear atom interfer- 
ometers with high brightness sources. 

VI. SUMMARY AND OUTLOOK 

Employing a spinor rather than a scalar condensate 
gives the opportunity to externally manipulate bright 
atomic solitons by conceptually simple magneto-optical 
methods without losing the stable solitonic behavior. 
Specifically, the characteristic properties of the gap soli- 
ton solutions lead to realistic manipulation schemes that 
have been demonstrated explicitly in a one-dimensional 
situation. This scheme has proven successful in exciting 
solitons with different velocities at a minimal atom loss 
rate as well as splitting, and reversing their direction of 
propagation. The combination of these techniques re- 
sulted in the demonstration of an atomic Mach-Zehnder 



interferometer, which might be of interest in atom optical 
sensors. We remark that although this was not explicitly 
discussed here, we have also found numerical evidence 
for stable three-dimensional gap solitons, when adding 
an optical potential for confinement in the transverse di- 
rection. 

It is now becoming amply evident that much of the fu- 
ture of atom optics lies in integrated systems, or " atom 
optics on a chip" iQ. Optical and magnetic waveguides 
and beam splitters have recently been demonstrated, but 
the coupling of a condensate into these guides remains a 
major experimental challenge. The use of bright solitons, 
with their potential to transport bright matter waves in 
a controlled fashion, might offer a solution to this prob- 
lem, and merits further investigations. Further topics of 
interest include the manipulation and control of bright 
solitons using phase imprinting methods that employ the 
AC~Stark effect, as have already been realized for dark 
solitons. 
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APPENDIX A: FOUR-PHOTON COUPLING 

As mentioned in the Introduction, the Raman coupling 
between the \Fg = l,mg = —1) and \Fg — l,mg = 1) 
Zeeman sublevels vanishes in the case of far-off resonance 
light, a result of destructive interference with other 3P3/2 
hyperfine excited states ||l^ . In this limit the system re- 
duces to an effective | = 1/2) ^ | Jg = 3/2) transition 
where no Am = ±2 transitions are allowed. This dif- 
ficulty can be overcome by adding a 7r-polarized laser 
beam from the side (e.g. along the X-axis), which will 
not transfer any momentum along the transverse direc- 
tion. A proper choice of detunings allows one to use 
another hyperfine ground state to mediate the Am = ±2 
transition, without populating it significantly. 

One disadvantage of this approach is that since this 
is now a four-photon process, the intensity of the lasers 
must be increased while still avoiding spontaneous emis- 
sion. Choosing a frequency difference between the tt- 
light and the cr-light a few MHz larger than the split- 
ting between the hyperfine ground states as shown in 
Fig. |ri|, three time scales govern the dynamics of the 
system. They are determined by (a) the coupling of the 
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ground state to the excited states, with detunings in the 
THz range; (b) the effective couphng between the hy- 
perfine ground states, with detunings in the GHz resp. 
MHz range; (c) the desired effective Am = ±2 cou- 
phng. As it turns out, couphng the \Fg = 2,mg = — f) 
and \Fg — 2, nig = 1) levels is a better choice than staying 
in the Fg — 1 manifold, as this minimizes the coupling to 
the outer states \Fg = 2, nig = ±2). Employing the inter- 
mediate \Fg — l,mg — 0) state makes the system aware 
of the hyperfine structure and leads to a non-vanishing 
coupling. 

Separating the time scales leads to an effective coherent 
evolution between the \Fg = 2, nig = ±1) states. It takes 
the same form as in Eq. (|]), but with 
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9eff 



1/192, 



(Al) 



where £^ is the amplitude for the cr-light resp. St^ for the 
TT-light and the detunings S and A are defined in Fig. 

The value of geff involves Clebsch-Gordan and 
Wigner-6J coefficients for all hyperfine transitions and 
is clearly a four-photon term due to the product of 
four electric field envelopes. Choosing the wavelength 
of the (T-light Xs = 985 nm as before, A = — Wo- = 
1.85 GHz and laser intensities 1^ — 100 KW/cm^ and 
= 158.3 KW/cm^ we obtain an effective coupling 
l/{geffS'^g) — 685.0 /is, as in Sec. ||. 

Numerical simulations of the full level scheme of 
Sodium is in exc ellent agreement with the derived cou- 
pling in Eq. (Al) and shows that the intermediate level 
\Fg — l,r7ig — 0) only contains around 0.1 percent 
of the atomic population. All other levels are negligi- 
bly populated, they are indeed far-off resonance. Since 
\Fg = 1, rrig = 0) is a ground state, spontaneous emission 
is not an issue. 

The intensities might be reduced if d were chosen 
smaller, but that would of course alter the properties 
of the soliton, since they depend on the wavelength of 
the optical fields. Note that we cannot apply permanent 
magnetic fields to remove the degeneracy of the hyperfine 
ground states and thus push non-participating states out 
of resonance since this would destroy the solitons. The 
light shifts induced by the lasers are the same for each 
manifold and will not remove this degeneracy. 
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FIG. 1. (7+ (T_ configuration: Counterpropagating laser 
fields induce far-oflf resonance Raman coupling between hy- 
perfine states. 
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FIG. 2. Total Number of atoms as a function of parameter 

rj for a stationary soliton and a rnoA-ing solitou with v — 0.83 
(i.e 1.47 cm/s). Moving solitons always contain less atoms 
than the stationary solitons. The transverse width was chosen 
to be 24.2 fxm, corresponding to an effective transverse area 
At = 919 nia^. 




FIG. 3. Width of the soliton as a function of parameter rj 
for a stationary soliton and a moving soliton with v = 0.83 
(i.e. 1.47 cm/s). For the regime of small r? we obtain widths 
of 40-100 urn. 
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FIG. 4. The propagation of a soliton excited at w = over 
around 50 soliton periods with initially 55000 atoms (shown is 
the total density) . It displays oscillations in the peak density 
due to imperfect initial conditions. Further numerical investi- 
gation shows that these oscillations will eventually damp out 
after around 1 s, accompanied by a loss of only two percent 
of the atoms from the soliton. 
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FIG. 5. The propagation of a moving soliton excited at a 
velocity of 0.18 cm/s over around 50 soliton periods (shown 
is the total density). The actual velocity seen in the plot is 
very close to the excitation value. Peak density oscillation are 
again due to imperfect initial conditions. 
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FIG. 6. The splitting of a condensate into two moving 
solitons, traveling into opposite directions (shown is the to- 
tal density). Their velocity is approximately ±0.9 cm/s, the 
propagation covers around 25 soliton periods. 



FIG. 9. Demonstration of an atomic Mach-Zehnder inter- 
ferometer: The initial condensate is split into two counter- 
propagating fast solitons, then their direction is reversed and 
they collide (shown is the total density). 
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FIG. 7. Two solitons 82.0 ms after the splitting. The plot 
reveals that each soliton has still a small fraction of the other 
state bound to it. 
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FIG. 10. Interference pattern when the two solitons collide 
in the Mach-Zehnder configuration: The total density (solid 
line) is symmetric, whereas the interference pattern of each of 
the two orthogonal Zeeman states (dotted and dashed lines) 
is asymmetric due to the shape of the colliding wavepackets. 
The contrast in the total density pattern is around 30%. 
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FIG. 8. The propagation of a fast soliton with velocity 
u ~ 1 cm/s whose direction is reversed after 100 ms. The 
process of reversing causes small loss: a packet travels on to 
the right. 
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FIG. 11. Four-photon configuration: Tlie \Fg = l,m = 0) 
level mediates the coupling between the \Fg = 2,m = ±1} lev- 
els. The red detuning of the 7r-light ensures that all other 
levels are off-resonance. Although the state \Fg = l,m = 0) 
is slightly populated, this ensures that we are sensitive to the 
hyperfine structure. 
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